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We discuss the possibility of coexistence of spin-singlet superconductivity and ferromagnetism in a 
model where the same electrons are assumed responsible for both of them. Our calculations include 
both zero and finite momentum pairing states with both s-wave and d-wave pairing symmetry. Under 
the mean-field approximation, the thermodynamic potential of the non-magnetic superconducting 
(SC) state is shown to be always lower than that of the superconducting ferromagnetic (SF) state. 
It follows that the spin-singlet SF state is energetically unfavorable, and a spin-triplet SF state is 
more likely to survive in metals such as UGe2 and ZrZn2. 
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I. INTRODUCTION 



Recently, the observation of superconductivity in fer- 
romagnetic metals, UGe2Q, ZrZnJ] and URhGoJ has re- 
newed the interest on the coexistence of ferromagnetism 
and superconductivity. The investigation on the su- 
perconducting ferraomagnets traces .back to the original 
wojjks of ClogstonQ, ChandrasekharEl, and Abrikosov et 
al.u. At the same time, the possibility of a finite momen- 
tum pairing state coexisting with the long range ferret 
magnetic order was also revealed by Fulde and FerrellEj 
and by Larkin and OvchinnikovO. This finite momentum 
pairing state|-iS|-usually called the FFLO state. These 
early worksffltffl'H focused on the superconductivity in 
the metals with a spin-exchange field, such as produced 
by ferromagnetically aligned impurities. In such a su- 
perconducting fcrromagnet, there are two kinds of elec- 
trons, respectively, responsible for ferromagnetism and 
superconductivity. One is the localized electrons forming 
a ferromagnetic background in the metal by way of an 
indirect exchange coupling through itinerant electrons, 
the other is the itinerant electrons forming the Cooper 
pair due to an effective attractive interaction. There- 
fore, the magnetic exchange energy / is independent of 
the superconducting gap A. If the spin-exchange field is 
weak enough, the superconductivity can appear against 
the ferromagnetic background. The ground state of the 
system is determined by the ratio 1/A.q, where Ao is 
the gap in a non-magnetic superconductor. In the three- 
dimensional s-wave case, for I/Aq < 0.707 all of the 
itinerant electrons near the Fermi level form the Cooper 
pairs with opposite spins and the center-of-mass momen- 
tum equal to zeroo, while for 0.707 < //A < 0.754 part 
of the itinerant electrons near the Fermi level form the 
Cooper pairs with a finite center-of-mass momentum and 
the unpaired electrons show a finite magnetization which 
is the paramagnetic response of the itinerant electrons 
to the exchange field caused by localized spinaJ'cl. For 
J/Ao > 0.754, the formation of the Cooper pairs is to- 
tally suppressed and the superconductivity is destroyed. 
In the two-dimensional d x i^ y i superconductor, with the 



increasing of the exchange field, the ground state of the 
system is changed from a zero momentum pairing state 
first to a finite momentum pairing state and then to a 
normal state, the same as in the s-wave case except that 
the zero momentum pairing d x ^. y 2 state in the magnetic 
field has a finite magnetization!! 

In the recently discovered superconducting ferromag- 
nets UGe2cl and ZrZnJa, the superconductivity and the 
ferromagnetism disappear at the same critical value, p c , 
under the application of the hydrostatic pressure. This 
feature may suggest that the same electrons are respon- 
sible for both ferromagnetism and superconductivity in 
these novel materials, in contrast to the conventional case 
of a metal with magnetic impunities. Very recently, along 
this direction, Karchev et altB have developed a theoret- 
ical model, in which the long range ferromagnetic or- 
der does not result from an indirect exchange coupling 
between the localized spins, but is a consequence of a 
spontaneously broken spin rotation symmetry of those 
itinerant electrons which participate in the Cooper pair 
formation. 

In this paper, we employ the model in Ref. 10 to dis- 
cuss the possibility of the coexistence of spin-singlet su- 
perconductivity and itinerant ferromagnetism. There are 
two kinds of interactions between the itinerant electrons 
in this model. One is an attractive interaction in the 
Bardeen-Cooper-Schrieffer (BCS) form resulting in the 
superconductivity, the other is an exchange coupling re- 
sulting in the ferromagnetic order. Thus the magnetic 
exchange energy and the superconducting gap are related 
to each other and will be solved self-consistently. Under 
the mean-field approximation, the thermodynamic po- 
tential of the superconducting ferromagnetic (SF) state 
and that of the non-magnetic superconducting (SC) state 
are obtained analytically at zero temperature. Our cal- 
culations include both the s-wave and d-wave cases. For 
each pairing symmetry, both the zero momentum and the 
finite momentum pairing states are discussed. It is shown 
that the thermodynamic potential of the SF state is al- 
ways higher than that of the non-magnetic SC state in all 
four cases. Therefore, the coexistence of ferromagnetism 
and spin-singlet superconductivity can not be realized if 
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the same electrons are assumed responsible for both of 
them. A spin-triplet superconductivity is more likely to 
survive in-those novel superpspducting ferromagnets such 
as UGezElBy and ZrZn 2 BEMl3. In Sec. II to Sec. IV 
we will explicitly illustrate this viewpoint for each case 
and in Sec. V our results are briefly summarized. 



II. ZERO MOMENTUM PAIRING s-WAVE 
CASE 

Our starting point is the model Hamiltoniar0 
H = J d*r £ c W_^V 2 -/,Wl 

- g J d 3 rc\(r)cl(r) Cl (r)ci(r), (1) 

where c a (f) are the fermion fields with spin a =t,4 , 
S = i ^2 a , c^Taa'Ca' is the spin field, r is the Pauli 
matrices, m* is the effective mass of the electron, and 
/i is the chemical potential. Hamiltonian ([j]) contains a 
kinetic energy term, a ferromagnetic exchange coupling 
term with strength J, and a four-fermion attractive in- 
teraction term with strength g. The coupling constant g 
has a finite value only in the thin shell of the width 2e c 
around the Fermi surface, as in the standard BCS theory. 
Under the -mean-field approximation, Hamiltonian (|l|) is 
reduced toEl 



Hmf = E e P (4t c pT + c li c pi 

E (4,^-4^0 + \jm 2 



p 

JM 
~2 



IAI 



(2) 



and 



M = ^E((4.i c p-i)-(4,T c ^T>) 5 (3) 



(4) 



where e p = p /(2m 



H is the band energy measured 
from the chemical potential, (• • • ) represents the thermo- 
dynamic average, M defines the magnetization of the sys- 
tem, and A is the superconducting gap. We note that the 
two constant terms in Eq. ^ result from the mean-field 
approximation, A 2 /g comes from the BCS interaction 
and JM 2 /2 from the exchange coupling. Here the mag- 
netization defined in Eq. (^) arises from a spontaneously 



breaking of spin rotation symmetry of the itinerant elec- 
trons, which is different from a paramagnetic response 
to a magnetic field caused by localized spins. There- 
fore, both the gap A and the magnetic exchange energy 
/ = JM/2 are determined by Eqs. (|) and (|) self- 
consistently, unlike in the conventional case of a metal 
with magnetic impurities, where th^ejSchauge energy is 
considered as an external parameterofflllLlla. 

By means of a Bogoliubov transformation, Hamilto- 
nian (^J) can be diagonalized as 

Hmf = E ( E p4 a P~ + tfPpPp) + E° (5) 
p 

with 

p 

where ap and [3p are the Bogoliubov fermion fields with 
excitation energies 



E f3 = Jel + \A\ 2 -I. 



(7a) 
(7b) 



The self-consistent equations (||) and (Q) take the form 

m = JE«-<)> ( g ) 



S|A| 



E 



(9) 



where n^' 13 is the momentum distribution function of the 
corresponding Bogoliubov fermions. At zero tempera- 
ture, rip' 13 is zero for ^ > and one for Ep"' 13 < 0. 
The thermodynamic potential at zero temperature can 
be obtained by averaging out the mean-field Hamiltonian 
(||). We take the thermodynamic potential of the ideal 
Fermi gas as the origin of the energy in the following 
discussions. 

Obviously, the self-consistent Eqs. (||) and @ have 
a non-magnetic SC solution with a finite gap Ao and 
zero magnetization. By replacing the summation over 
the momentum p in Eq. (|j) with an integral over e p , for 
the weak coupling limit gN(0) <C 1, one can easily find 
that the energy gap in the non-magnetic SC solution is 
given by 



Aq = 2e c exp 



gN(0) 



(10) 



where -/V(O) is the density of states (DOS) at the Fermi 
level. Therefore, the thermodynamic potential of the 
non-magnetic SC state at zero temperature takes the 
well-known form 



sc 



iiV(0)A 



(11) 
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Next, we discuss the SF solution at zero temperature 
with both finite gap and finite magnetization. Here, the 
gap is a real number, and the self-consistent Eqs. (||) and 
M turn to 



M 



N(0) I de, 



(12) 



for J/Ao < 0.707 the itinerant electrons form a non- 
magnetic SC state in the presence of an exchange field of 
the ferromagnetic background. In the present case, the 
external exchange field is absent, so is the spontaneous 
magnetization of the itinerant electrons. As a result, the 
coexistence of s-wave superconductivity with zero mo- 
mentum pairing and itinerant ferromagnetism can not 
be realized. 



gN(0) 



In 



2ec 
A 



de 

V/^-a^ Ve 2 + A 2 



(13) 



The integrals over e p in the self-consistent equations are 
confined in a very thin shell around the Fermi level, and 
therefore we have replaced the DOS in the integrals by 
iV(0). It is shown that, in the SF state, only part of 
the itinerant electrons whose energies are in the range 
\JP- — A 2 < |e| < e c form the Cooper pairs while other 
electrons with energies in the range < |e| < \/ 1 2 ~ A 2 
remain unpaired. These unpaired itinerant electrons give 
rise to a finite spontaneous magnetization. Completing 
the integrals in Eqs. ( |l2|) and (|l3|), one finds that, for 
r > 1, the SF solution is given by 



A 



: A , 



(14) 



where r = JN(0)/2 is the dimensionless measurement 
for the exchange coupling strength. Part of the electrons 
near the Fermi level remain unpaired so that the super- 
conducting gap A in the SF state is always less than Ao. 
Here, the exchange energy / results from the exchange 
coupling between the itinerant electrons and is related 
to the gap A, in contrast to the case of a metal with 
impurities. For fixed r, the ratio I/A is a constant, in- 
dependent of A. When the gap vanishes the exchange 
energy also decreases to zero. Substituting the SF solu- 
tions for the magnetization and the gap into Eqs. (|^) 
and (JgJ) and averaging out the mean-field Hamiltonian 
(|J), one obtains the thermodynamic potential of the SF 
state at zero temperature as 



-iV(0)A 2 



(15) 



Comparing Qsc with Ogp, one finds that the SF state 
in the zero momentum pairing s-wave case is always en- 
ergetically unfavorable due to A < Ao- This result can 
be understood by the following argument. The sponta- 
neous magnetization in the SF state is very weak and 
the exchange energy / in the SF state is always less than 
Ao. From Eq. (0) and (|TJ) one finds that the Zeeman 
energy gained by the unpaired electrons can not com- 
pensate for the loss of the condensate energy in depair- 
ing the Cooper pairs. Therefore, all the itinerant elec- 
trons near the Fermi level remain paired with opposite 
spins, and the non-magnetic SC state is always more sta- 
ble than the SF state. This feature is similar to that in 
the conventional case of a metal with impurities, where 



III. FINITE MOMENTUM PAIRING s-WAVE 
CASE 

In the conventional case of a metal with impurities, the 
formation of a condensate of finite center-q£-|Hiass mo- 
mentum turns out to be more advantageousDB. Such a 
finite momentum pairing state, where part of the itiner- 
ant electrons near the Fermi level form the Cooper pairs 
with a finite center-of-mass momentum and the unpaired 
electrons show a finite magnetization, is the ground state 
of the system when the exchange energy / is in a proper 
range. However, if the exchange energy is caused by the 
spontaneous magnetization of the itinerant electrons and 
thus related to the gap, the conclusion is completely dif- 
ferent, as will be shown below. 

Here, the center-of-mass momentum of the Cooper 
pair, <f, is not equal to zero, and thus the supercon- 
ducting gap, turns out to be a periodic function of the 
coordinatescl, such as A(r) = Aexp (iq ■ r). The excita- 
tion energies of the Bogoliubov fermians in such finite 
momentum pairing state are given byc3 



E° = 



= 

p 



1 



p.q 



A 2 + 1 + -v F qcos6, (16a) 
1 



p.q 



A 2 - I- -v F qcos6, (16b) 



where e„ 



= (e p - +m + e_ f7+?72 )/2 = (p 2 + q 2 )/(2m*), 
vf is the Fermi velocity, and 9 is the angle between the 
momentum p and q. The diagonalizcd Hamiltonian, the 
self-consistent equation of the magnetization and the gap 
equation have the same forms as Eqs. (||), (||), (|J), and 
(^), respectively, provided that e„ is replaced by e Piq and 
the excitation energies in Eqs. (Ra) and (f?]b) are replaced 
by those in Eqs. (|l6|a) and (|16[b). The momentum q is 
determined by the minimization of the thermodynamic 
potential, and VFq/2 is at most of the order A o 0. By 
replacing the summation over the momentum p with an 
integral over e = e Pj q/A, the self-consistent equations at 
zero temperature can be written as 



hi- 



dxde 



An 



4 \Je°<o ^ 
-Fx{l',y), 



dxde 



1 Je%<o Ve 2 + 1 



(17) 




dxde 



ES<0 



dxde =rF 2 (l',y), 



(18) 
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1 



on 



N(0)A 2 dy 



F 3 (l',y) = 0, 



(19) 



where I' = I/A, y = (vpq/2)/ A, x = cos9 and -Fj.,2 3 are 
functions of /' and y, independent of A. From Eqs. fll8| ) 
and (|l9| ) , one finds that J' and y are determined only by 
r. Hence, the ratio A/Ao in Eq. (|l7| ) is also a constant 
for fixed r. Combining with Eq. ([Uj), one obtains 



51) 
dA 



-iV(0) 



1 



(20) 



This result is quite different from that in the conventional 
case where / is an external parameter and the right-hand 
side of Eq. (|l7|), for fixed /, depends on the gap explic- 
itly. The thermodynamic potential of the SF .state with 
moving pairs is given by applying the formula!!!] 



n SF - n(A = 0) = 



A ; 



:__9_ 

dA 



dA, 



(21) 



where f2(A = 0) is the energy when the gap vanishes. 
In our cases, the ratio I' = I/A is a constant for fixed 
r. When the gap vanishes the exchange energy also de- 
creases to zero. Therefore, f2(A = 0) here is actually the 
energy of the ideal Fermi gas and is set to zero in our 
discussions. Completing the integral in Eq. (pl[), one 
obtains 



fi SF = -^V(0)A 2 . 



(22) 



The energy in the SF state with moving pairs defined 
in Eq. ( f22[ ) takes the same form as that with immobile 
pairs defined in Eq. (]l5|) , but the gap functions in the two 
cases are different, because they are determined by two 
different self-consistent equations. Noting that the gap A 
is always less than Ao, we find that the thermodynamic 
potential of the SF state with moving pairs is still higher 
than that of the non-magnetic SC state. Our result is 
different from that in FFLO state, in which part of the 
itinerant electrons near the Fermi level remain unpaired 
and give rise to a finite magnetization in an exchange field 
of 0.707 < I/Aq < 0.754 and other electrons form the 
Cooper pairs with a finite center-of-mass momentum!]. 
This is because that the Zeeman energy gained by the 
unpaired electrons in the SF state in the present model 
is less than that in FFLO state. In the conventional case, 
the exchange field is caused by localized spins and the ex- 
change energy I is an external parameter. Each itinerant 
electron moving in the exchange field has an additional 
energy / or —I depending on the spin orientation of the 
electron. The difference in number between spin-up and 
spin-down electrons results in a finite magnetization M. 
The Zeeman energy of the system is given by —2IM. In 
the present model, the exchange field is due to the spon- 
taneous magnetization of the itinerant electrons and the 
exchange energy / is determined self-consistently. There- 
fore, the additional energy / gained by the electron in the 
exchange field here is actually the interaction energy be- 
tween this electron and other itinerant electrons. The 



Zeeman energy of the system is obtained by adding the 
exchange energy of every electron together, but the inter- 
action energy between the two electrons should be coun- 
tered only once. Thus, the Zeeman energy here is —IM, 
which is only a half of that in an external field with the 
same /. This feature is embodied in Hamiltonian (Q) via 
the constant term JM 2 /2. Such a loss in the Zeeman en- 
ergy results that the SF state with moving pairs is also 
energetically unfavorable and the coexistence of the itin- 
erant ferromagnetism and the superconductivity is more 
difficult. As a result, if the exchange energy is caused 
by the spontaneous magnetization of the itinerant elec- 
trons and determined self-consistently, the s-wave non- 
magnetic SC state is always more stable than the s-wave 
SF state, no matter whether the center-of-mass momen- 
tum of the Cooper pair is zero or not. 



IV. d-WAVE CASE 

A. Zero momentum pairing d-wave case 

We take a two-dimensional d x 2_ y2 superconductor for 
example. The gap function takes the form A(9) = 
Ad cos(2#), where 9 is the azimuthal angle of momentum 
p. The self-consistent equation for the magnetization re- 
mains the form of Eq. (M) and the gap equation is given 
byl 



E 



e 2 + A 2 , cos 2 (29) 



:Cos 2 (26>), (23) 



where and are the momentum distribution func- 
tions with the energies 



E; = ^e 2 +A 2 cos 2 (2#)+/, (24a) 
El = Jd + Al cos 2 {29) ~ I. (24b) 



At first, we consider the non-magnetic SC solution 
where 1 = and thus = = 0. Carrying out the 
summation in Eq. (]23|), the gap in the non-magnetic SC 
state is obtained by 



A d0 = 2.43e c exp(- 



gN(0)' 



(25) 



By using Eq. (pl|), the thermodynamic potential of the 
_. y 2 non-magnetic SC state at zero temperature is 



given by 



ttsc = - i ^V(0)A 2 o . 



(26) 



Next, we turn to the SF solution, where only part of 
the itinerant electrons near the Fermi surface form the 
Cooper pairs while other electrons whose energies sat- 
isfy the inequality, Eg < 0, remain unpaired. Replacing 
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the summations over momentum p in Eqs. (||) and ( pjj ) 
by the integrals over the energy e = e p /Ad and the az- 
imuthal angle 9, the gap equation and the self-consistent 
equation of the magnetization are reduced to 



In 



Vio 



1 



= / d9de 



cos 2 (26») 



-ES<0 



\/e 2 + cos 2 (26>) 



(27) 



Combined with Eqs. ( |2q ) and (|21|), the thermodynamic 
potential of the SF state with finite momentum pairing 
is obtained as the same form as defined by Eq. ( p9| ) but 
the gap here is different from that in the zero momentum 
pairing case. It then follows that the energy in the d x 2_ y2 
SF state with moving pairs is also higher than that of the 
non-magnetic SC state due to Ad < A d o and therefore 
the SF state can not be realized. 



I 

Ad" 



47T 



dOde. 



E%<0 



(28) 



From Eq. (28), it follows that the ratio //Ad is a con- 
stant for fixed r, and therefore the ratio Ad/Ado in Eq. 
(^7]) is also a constant for fixed r. The exchange energy 
I is related to the gap Ad, and decreases to zero as the 
gap vanishes. With the help of Eqs. (21) and (p5|), the 
thermodynamic potential of SF state is given by 



(29) 



The gap in the SF state is always less than that in the 
non-magnetic SC state. Hence, the d x 2_ y2 SF state al- 
ways has higher energy than the d x 2_ y2 non-magnetic SC 
state. This result is the same as that in the s-wave case. 



B. Finite momentum pairing d-wave case 



In the d x 2_ y2 SF state with finite momentum pair- 
ing, the excitation energies of the quasi-particles take the 

form E% = 
and E$ = w t p q 



yj^ q + Ajj cos 2 (20) +I+(v F q/2) cos(0 - 9 q ), 



A 2 cos 2 (26>) - 1 - (v F q/2) cos(6> - 9 q ), 
where 9 q is the azimuthal angle of the momentum q which 
is determined by the minimization of the thermodynamic 
potential. The self-consistent equations in the finite mo- 
mentum pairing d-wave case have the same structure as 
that in the s-wave case. The equations for the mag- 
netization and the momentum of the pair can be re- 
duced to three expressions which only depend on //Ad, 
(i>F<z/2)/Ad, 9 q and r. Therefore, the variables //Ad, 
(vFq/2)/Aa and 9 q are found to be constant for fixed r. 
The gap equation is reduced to the expression which only 
contains the variables Ad/Ado, //Ad, (vpq/2)/ 'A<i, and 
9 q . Thus, the ratio Ad/Ado is also a constant for fixed r. 



V. SUMMARY 

In this paper, we employ a model including both an 
attractive interaction and an exchange coupling between 
the itinerant electrons to study the possibility of the co- 
existence of spin-singlet superconductivity and itinerant 
ferromagnetism. In this model the electrons responsible 
for the ferromagnetism and those forming the Cooper 
pairs are the same and thus the exchange energy between 
the two spin subbands is determined self-consistently, as 
opposed to the conventional case where the exchange en- 
ergy is an external parameter. Under the mean-field ap- 
proximation, the self-consistent equations of both the su- 
perconducting gap and the exchange energy are consid- 
ered simultaneously, and the thermodynamic potential in 
the non-magnetic SC solution and that in the SF solution 
are obtained analytically at zero temperature. We dis- 
cussed four cases including both the zero and finite mo- 
mentum pairing state with both the s-wave and <i-wave 
pairing symmetry. One finds that the exchange energy 
and the superconducting gap in the SF solution are re- 
lated to each other. It is shown that the Zeeman energy 
gained by the unpaired electrons in the SF state can not 
compensate for the loss of the condensate energy in de- 
pairing the Cooper pairs and thus the thermodynamic 
potential of the SF state is always higher than that of 
the non-magnetic SC state in both the s-wave case and 
the d-wave case, no matter whether center-of-mass mo- 
mentum of the pair is zero or not. Therefore, the coexis- 
tence of ferromagnetism and spin-singlet superconductiv- 
ity can not be realized if the same electrons are assumed 
responsible for both of them. The present results indi- 
cate that a spin-triplet superconductivity is more likely 
to survive in those novel superconducting ferromagnets 
such as UGe2-apd-Zi^n2, which are consistent with some 
new progress 1 1 
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